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Abstract 

The canonical quantization of homogeneous cosmologies is considered in the 
high anisotropic limit. Exact wavefunctions are found in this limit when the 
momentum constraints are reduced at the classical level. Lorentzian solutions 
that represent tunnelling from classically forbidden regimes are identified. So- 
lutions to the modified Wheeler-DeWitt equation are also found for the vacuum 
Bianchi IX model when a quantum reduction of the momentum constraints is 
considered. 
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Quantum cosmology applies the concepts of quantum mechanics to the Universe as 
a whole. The majority of studies in this field invoke the Dirac quantization procedure 
[IJ. In this approach the wavefunction of the Universe is annihilated by the operator 
versions of the Hamiltonian and momentum constraints || . However, it is not known 
how to solve these constraints in full generality. Indeed, the configuration space 
of the wavefunction is infinite-dimensional. In practice, therefore, one invokes the 
'minisuperspace' approximation and applies the procedure only to those cosmologies 
that represent homogeneous solutions to Einstein's field equations @, [|. These are 
the Bianchi models and the Kantowski-Sachs Universe. Hence, the problem is reduced 
to a quantum mechanical system with a finite number of degrees of freedom. 

It is not clear whether the results derived from this minisuperspace quantization 
represent a valid approximation to the full theory of quantum gravity 0, |5j . However, 
the main justification for this approach is that it provides a solvable framework in 
which the problems associated with quantum cosmology may be addressed. For ex- 
ample, there are difficulties in extracting physical predictions from the wavefunction 
and there are also problems with defining a Hilbert space of states. (For a review see, 
e.g., Refi). 

One may adopt one of two approaches when performing a canonical quantization 
of homogeneous Universes. The simplest method is to reduce the momentum con- 
straints at the classical level before quantization and most investigations to date have 
adopted this view. The Wheeler-DeWitt equation represents the quantum Hamil- 
tonian constraint and governs the evolution of the wavefunction on minisuperspace 
0. However, a more rigorous and consistent method is to reduce the momentum 
constraints at the quantum level. In this way all constraints are treated in a similar 
fashion. These two approaches are not necessarily equivalent and can lead to different 
results 0, §. 

In either case, however, very few exact wavefunctions have been found to date. 
The purpose of this letter is to derive new families of solutions to the quantum 
constraints of the Bianchi class A models. We consider wavefunctions on a region of 
minisuperspace corresponding to the high anisotropic limit of these models. 

We first investigate the approach whereby the momentum constraints are reduced 
prior to quantization. We show that the resulting Wheeler-DeWitt equation may 
be written in the form of the unit-mass Klein-Gordon equation. We then solve the 
Wheeler-DeWitt equation for the vacuum Bianchi IX model when a full quantum 
reduction of the momentum constraints is performed. 

The homogeneous Bianchi space-times have a topology R x G3, where G3 is a 
three-dimensional Lie group of isometries transitive on space-like three-dimensional 
orbits §J. The world-interval of the space-time is 

ds 2 = -dt 2 + h ab u V, a, b =1,2, 3, (1) 

where the metric on the surfaces of homogeneity is a function of t alone and u a are 
one-forms. The isometry of the three-surface is determined by the structure constants 
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Type 


I 


II 


VI 


VIIo 


VIII 


IX 


m ab 


(0,0,0) 


(1,0,0) 


(1,-1,0) 


(1,1,0) 


(1,1,-1) 


(1,1,1) 



Table 1: The diagonal components of m for each Bianchi type in the class A. 



C a b c of the Lie algebra of G 3 . These may be decomposed as C a b c = m ad edbc + <$° [6 a 



c , 



where a c = C a ac and m ab is symmetric fllP]. The Jacobi identity C a bf c C b d e i = is 



satisfied if and only if a& is transverse to m , i.e. m a>b = 0. The Lie algebra belongs 



to the Bianchi class A if a& = and to the class B if ab ^ p0| . It is well known 
that the evolution of class B models cannot be described in terms of a standard 
Hamiltonian treatment. This is related to the fact that these models cannot admit a 
spatially compact topology 0. Consequently, we restrict our attention to the class 
A. This consists of Bianchi types I, II, VIo, VIIo, VIII and IX. The Lie algebra of 
each type is uniquely determined up to isomorphisms by the rank and signature of 
m ab . 

We will assume that the matter source is a single, massless, minimally coupled 
scalar field <p. The classical dynamics of these cosmologies is determined by the 
time-space and time-time components of the Einstein field equations. These may be 
expressed as the momentum and Hamiltonian constraints, respectively. The former 
is given by 

P d = rr a c m cb e dab = 0, (2) 

where the conjugate momentum variable 7i ab represents the integral of the momentum 
density over the spatial hypersurfaces. The Hamiltonian constraint, on the other 



hand, takes the form [11] 



rr ab 7r ab - ^tt 2 ) + V 2 (m ab m ab - ^m 2 ) + ^vr 2 = 0, (3) 

where 7T^ is the momentum conjugate to the scalar field, V represents the volume of 
space and indices are raised and lowered with h ab and h a b, respectively. 

At the classical level, the momentum constraints (0) imply that -K ab , m ab and h a b 
may be simultaneously diagonalized on a given spatial hypersurface of constant t. It 
then follows directly from the field equations that these quantities are diagonal on all 
other hypersurfaces. Hence, one may diagonalize the variables at the classical level 
without loss of generality and the momentum constraints therefore become trivial. 

In this case one may choose a basis where each diagonal component of m ab is ±1 
or 0. The six cases that constitute the Bianchi class A are shown in Table 1. The 
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three-metric is then written as 



hab(t) 



„2a(t) 



ab 



(4) 



where the matrix (3 ab = diag (3 + + a/3/3_, /3 + — \/3/3_, — 2/3 + represents the degree of 
anisotropy in these models. The parameter e a may be viewed as an averaged scale 
factor of the Universe. The Hamiltonian constraint (|3|) then reduces to 



— 7T„ 



7T 



■Kr. 



+ 7ci + u = o 



where iij (j = a,(3±,<f>) are the conjugate momenta, 



(5) 



U= 3 



(mn/iu) 2 + {m 22 h 22 f + (m 33 h 33 Y 



-2m xx m 22 h xx h 22 - 2m 11 m 33 h l ih 33 - 2m 22 m 33 h 22 h 



(6) 



represents the 'superpotential' and m ab = diag [mu, m 22 , m 33 ] p. (The direct depen- 
dence of the superpotential on V is eliminated by performing a linear translation on 
a). The superpotential for the Bianchi type I cosmology vanishes identically and we 
shall not consider this model further. 



Quantization of these cosmologies follows by imposing the algebra [j, ttj 



i. 



The Hamiltonian constraint (^) is then promoted to an operator that annihilates 
the state vector \1/ of the Universe. The result is the Wheeler-DeWitt equation on 
minisuperspace: 



d 2 



d 2 



d 2 



da 2 d(3l d(3 2 _ 



d 

da 
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+ u 



^ = 0, 



(7) 



where the constant p accounts for ambiguities that arise in the operator ordering [[12 

It is convenient to rescale the wavefunction by \P = \^e iq ^~ pa , where q is an 
arbitrary, complex constant. Eq. (^) then simplifies to 



D 2 



d 2 



d 2 



da 2 dm d(3l 



V 



u 



m = o. 



A unified exact solution to this equation can be found for all Bianchi types if we 
choose the factor ordering p 2 = 9 + q 2 . It can be verified by direct substitution that 
the wavefunction has the form 



where 



v|> = e (3-p)a+iq<f> e -i 



±7:rn ab h ab 
o 



(9) 



(10) 



and summation over indices is implied. This solution generalizes the exact solution 
found previously for the type IX model ||13||. In the vacuum case (q = 0), the function / 
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is itself a solution to the Euclidean Hamilton- Jacobi equation. It therefore represents 
a Euclidean action for these cosmologies. 

In general, an oscillating solution to the Wheeler-DeWitt equation may be inter- 
preted as a classically allowed, Lorentzian geometry ||I4|. A cosmological singularity 
then arises when the wavef unction undergoes an infinite number of oscillations. On 
the other hand, a classically forbidden Euclidean geometry corresponds to a non- 
oscillatory, exponential wavefunction. In particular, if the wavefunction is suitably 
damped at large three-geometries and remains regular when the spatial metric degen- 
erates, it may be interpreted as a quantum wormhole [I5| . Solution (|D|) is an example 
of this latter type of wavefunction and therefore represents an anisotropic quantum 
wormhole. 

Further unified solutions to Eq. @ can be found in the limit where f3± ^> 1. 
This region of minisuperspace is interesting because quantum fluctuations of the 
gravitational field in the early Universe may have introduced considerable anisotropy 
into the initial conditions. However, observations indicate that the present Universe 
is highly isotropic. It is possible, therefore, that a better understanding of this high 
anisotropic regime may provide insight into how the Universe evolved into its current 
state. 

It follows that ^22/^11 *C 1 and h^/fin <C 1 when /3± > 1. The first term in Eq. 
(§) therefore dominates the superpotential. Formally, this is equivalent to choosing 
mn = 1 and 17122 = 1^33 = 0. Hence, this particular limit of the Bianchi class A may 
be investigated by considering the type II model. The Wheeler-DeWitt equation (§) 
simplifies to 
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d 2 



d 2 



da 2 d(3 2 d(3 2 



+ q 2 



p 2 _|_ I e 4a+4/3 + +4 v / 3/3- 
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* = 



and may be solved by defining a new wavefunction 

$(a,Afc) = ^(a,/?±)e- r(v/5/3+ - /3 - )/2 , 
where r 2 = q 2 — p 2 . Substitution of this ansatz into Eq. ([11]) implies that 



(12) 



d 2 



2 



d 2 



da 2 dp 2 dp 2 



J^r-?— + r-^— + Ip4a+4/3 + +4 v / 3/3- 



$ = 0. 



(13) 



We now perform a change of variables to {u, v, w}, where 



U =l2 e 



v = — e 
12 



6a+3/3 + +3V3f3- 



-2a+f3++V3/3- 



(14) 



and w = w(a,j3±) is an arbitrary, twice continuously different iable function of the 
minisuperspace coordinates. We then search for wavefunctions that are independent 
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of this variable. In this case, Eq. 
Gordon equation 



(|13|) simplifies to the canonical, unit-mass Klein- 
d 2 



dudv 



$ = 0. 



(15) 



Hence, the {u, v} variables may be viewed as null coordinates over a region of 
(1 + l)-dimensional Minkowski space-time spanned by the space-like coordinate X = 
u + v and time-like coordinate T = u — v. These null coordinates are restricted to lie 
in the range (u, v) G (0, +oo), so X > \T\. We may therefore view the wavefunction 
$ as a classical particle of unit mass that is confined within the Rindler wedge of 
(1 + l)-dimensional Minkowski space-time. 

The semi-classical limit of the vacuum models may also be analyzed by introducing 
the variables fll4l). In the WKB approximation, one substitutes solutions of the form 
$wkb ~ e~ lS l n into the Wheeler-DeWitt equation and considers the limit h — > 0. 
Applying this substitution in Eq. (p~5|) implies that the Hamilton- Jacobi equation 
takes the form S U S V = —1, where a subscript denotes partial differentiation. One 
solution to this equation is S = —i(cu + c _1 t>), where c is an arbitrary, complex 
constant. However, the family of wavefunctions $ c = e~ cu ~ v ^ c are also exact solutions 
to the full Wheeler-DeWitt equation (|T5| ) [filj. Hence, the WKB approximation is 
exact in this case. 

|$ c | is bounded everywhere when Re c > and the wavefunction is exponentially 
damped for large a if Im c = 0. Unfortunately, however, it does not satisfy the quan- 
tum wormhole boundary conditions because it decays too rapidly ||15|| . On the other 
hand, more general solutions to Eq. (|15"D may be generated in terms of linear super- 
positions of this family. The class of wavefunctions $ c may therefore be physically 
relevant. In general, the superpositions have the form 



c 



dcM(c)e 



-cu—v/c 



(16) 



where M(c) is an arbitrary function of the parameter c and C represents the contour 
of integration in the complex plane ||16|| . 

Different solutions to Eq. ([15]) correspond to different choices for M(c) and C. 
For example, if we specify M(c) = ^ te ~ 3 ^ 3 , where e is an arbitrary, real constant, 
and perform the integration over the positive half of the real axis, we find that 
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dcc {i€ - 3)/3 e- cu - v/c 



K, 



ie/3 



(2vW) 



K 



ie/3 



l e 2a+2/3 + +2 v / 3/3- 

6 



exp 



~ ( 8a + 2(3+ + 2V3/L 
o v 



(17) 



where K is the modified Bessel function of order ie/3 [IT 
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The argument of this function is given by 



2\/uv 



6 6 



+2/3+ +2^/3- 



and corresponds to the high anisotropic limit of the Euclidean action ([101) . It is also 
directly proportional to the square root of the superpotential (§]). The modified Bessel 
function is exponentially damped for sufficiently large arguments and the wavefunc- 
tion takes the form $ e oc e - ' 7 ' for y/uv > |e|/6. In this region of minisuperspace, 
therefore, the wavefunction reduces to the Euclidean form of Eq. ( [Top . 

The nature of solution ([17]) is different for smaller values of the scale factor. Indeed, 
it oscillates for < |/| < \e\/3, so the boundary I = |e|/3 represents a point of 
maximum expansion. The wavefunction has the asymptotic form $ e oc u~ lt ^ 3 as 
a — > — oo and therefore represents plane waves in the minisuperspace variables when 
e 7^ 0. It oscillates an infinite number of times as the spatial metric degenerates and 
this behaviour is interpreted as a cosmological singularity. 



Solution ([TBI) may provide insight into the quantum nature of cosmological sin- 
gularities in these highly anisotropic Universes. The singularity arose in this case 
because a particular superposition of Euclidean solutions was considered. In view 
of this, it is natural to investigate whether the singularity may be eliminated by 
considering alternative superpositions. 

In order to pursue this possibility further, we introduce the new coordinate pair 



v 



This change of variables transforms Eq. (fL5[) into 



v 

v = — 

2 



2u. 



& 2 



cf 2 



$ = 



<9/i 2 dv 2 

and the general, separable solution to this equation is given by 

$ m = [ciAi(m + ji) + c 2 Bi(m + //)] [c 3 Ai(m + v) + c 4 Bi(m + v)\ , 



(19) 



(20) 



(21) 



where Ai(x) and Bi(a;) are Airy functions and {m,Cj} are arbitrary constants. 

This solution may be expressed in the form of Eq. ( |16D when Cj satisfy appropriate 
conditions. If we specify 



M(c) 



V2i 



exp 



2 

3? 



2m 

c 



(22) 



Eq. (jl6|) takes the form 



dc 



c c 



1/2 



exp 



-^ + (/i + z/ + 2m)^ + l(/i-z/) 2 



(23) 
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where c = — 2/c. Halliwell and Louko have shown how integrals of this form may 



be evaluated [17]]. Different choices for the contour C result in different products of 
Airy functions. Since we are interested in superimposing bounded wavef unctions, we 
assume that Re c > 0. We therefore choose the contour of integration to lie to the 
left of the origin. The result of the integration is 

$ m = Ai(/i + m)M(u + m) (24) 



and, modulo a constant of proportionality, this is equivalent to Eq. with c 2 = 

c 4 = 0. 

This solution has an interesting feature. The Airy function Ai(x) is exponentially 
damped for large positive arguments, but exhibits oscillatory behaviour if x < ]|i~8 



However, the variable (/i + m) is positive-definite for all values of the scale factor 
if m > 0. In this case, it follows that the wavefunction will represent classically 
forbidden geometries when [y + m) > 0, but will correspond to Lorentzian solutions 
when (u+m) < 0. Now, when the spatial volume of the Universe becomes vanishingly 
small (a — > —00), u — > and v — > +00, so fi — > v — > +00. Hence, the wavefunction 
vanishes, but does not oscillate, as the spatial metric degenerates. Consequently, there 
is no singular behaviour at the origin. As the scale factor grows, however, [y + m) 
decreases and eventually becomes negative. At this point the wavefunction begins 
to exhibit oscillatory behaviour. In effect, the Universe tunnels from a classically 
forbidden regime into a Lorentzian domain when v « — m. 

The above solutions were derived after the momentum constraints had been re- 
duced at the classical level. However, it is more accurate to view these constraints as 
operators that annihilate the wavefunction of the Universe. Consequently, it is better 
not to impose the assumption of diagonality before quantization. Indeed, Ashtekar 
and Samuel have argued that this restriction may result in a loss of generality at 
the quantum level 0. In the remainder of this work, therefore, we shall consider the 
quantization procedure when the assumption of diagonality is dropped. In particular, 
we will investigate the vacuum Bianchi IX model. 

The quantum constraints for this model were recently derived by Higuchi and 
Wald ||. We briefly summarize their main results here. They parametrized each 
point in minisuperspace in terms of an orthogonal matrix O and the eigenvalues 
{Ai, A 2 , A 3 } of the metric h. These eigenvalues satisfy the conditions Ai > A 2 > A 3 > 
and A1A2A3 = 1. The matrix O rotates the eigenvectors of the metric into a new 
orthonormal basis. It can be shown that the quantum momentum constraints are 
satisfied if the wavefunction is invariant under the action of OhO -1 . It must therefore 
be symmetric under permutations in \. 

If one chooses Ai = e 2 ^++^-\ A 2 = e 2 ^+~^^ and A 3 = e~ 4f3 +, the Wheeler- 
DeWitt equation may be written as 



r\2, 1 O Q 

da 2 C((3±) fr± d/3j d/3j 



^ = 0, (25) 
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where 



sinh (2V3P-) sinh (3/5+ - V%P-) sinh (3/3+ + y/30_ 



(26) 



and £ is a numerical constant. The non-trivial contribution from C{(3±) arises be- 
cause the volume element on superspace is affected by a term that depends on the 
eigenvalues of the metric. 

In this representation, the momentum constraints are satisfied if the wavefunction 
is invariant under 120-degree rotations in the (/?+, (3 J) plane and reflections in f}_. An 
appropriate linear combination of wavef unctions satisfying these symmetry conditions 
can always be constructed. Formally, one may write the full wavefunction as \I/ = 
$ + + + $(/?_ — > —/?_), where $ represents a particular solution to Eq. ( pop 
and R is the 120-degree rotation matrix in the (/?+,/?_) plane [19]. On a practical 



level, therefore, one need only find solutions to the modified Wheeler-DeWitt equation 



However, it is not clear at present how one might proceed to solve this equation in 
full generality. In view of this, we shall consider the two limiting cases where (3± ^> 1 
and (3+ > 1, < 1. 

In the former case C{(3±) « e 6/3++2v ^ /3 - 
simplifies toQ 



and U = 48m>. Eq. (2^) therefore 



d 2 



d 2 



d 2 



da 2 dp\ d{3 2 _ 
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d 



dp, 



2V3 







dp. 



£/-90f 



* = 0. 



(27) 



A quantum wormhole solution to this equation may be found by rescaling the 
wavefunction such that (p = \&e 3l3++ ^ l3 ~ . This wavefunction satisfies an equation 
that is formally equivalent to Eq. (B), where q = and p 2 = 90£ — 12. Hence, one 



(2? 



solution is given by <p = e 2v ^ for the special case £ = 21/90. 

Further solutions may be found by defining a new wavefunction $: 



where 



$(a,/?±) = V(a,P±)e 



73 



-C+P++C-P- 



12 - 90£ - 3 



(29) 



The Wheeler-DeWitt equation 

d 2 d 2 d 2 



transforms to 



da 2 dp\ d(3 2 



- V3 (12-90£) 



1/2 



d_ 



+ (12 - 900 



1/2 







4a+4/3++4 v / 3/3- 



dp. 
$ = 



(30) 



■"■Note that we are assuming implicitly that /3+ > /3_ in this analysis. This condition ensures that 
C(/3±) nas an approximately exponential form. 
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and this equation is formally equivalent to Eq. ([13]) with r = (12— 90£) 1 ^ 2 - It therefore 
transforms into the unit-mass Klein-Gordon equation ([TJ) when the null variables 



(fj) are introduced. Hence, the above analysis and solutions will also be relevant in 
this more general quantization procedure if the model is sufficiently anisotropic. The 
range of j3± must be such that C(j3±) can be viewed as a purely exponential function. 

We will conclude by considering the case where P + 3> 1 and |/3_| -C 1. It fol- 
lows that C(f3±) « 4a/3|/3_ |e 6/3+ in this limit and the Wheeler-DeWitt equation (g5|) 
therefore takes the approximate form 



d 2 



d 2 



d 2 



1 d 



da 2 dpi dp 2 (3. dp. 



+ (9 - 900 + e 4a+4/3+ /3! 



6 = 0, 



(31) 



where B = e 3f3+ ^. Without loss of generality, we have eliminated the numerical 
constant in front of the superpotential via a linear shift in a. 
To proceed, we change variables to the null coordinates 



In this case, Eq. 



p = a + p+, 
transforms to 



4- 



d 2 



d 2 



1 d 



dpdri dpi p. dp. 



r] = a — P^ 



+ (9 - 900 + P-^ P 



9 = 



(32) 



(33) 



and it follows that the wavefunction G is an eigenstate of rj. In order to solve this 
equation, we assume that G has the generic form [|J 



G = exp -iEr] - B{p) - XA(p)p 2 _ 



(34) 



where A and B are arbitrary functions of the null variable p and {A, E} are arbitrary 
constants. After substitution of this ansatz, Eq. (|33| ) separates into two, first-order 
ordinary differential equations: 



AiEB' + AXA + 9 - 90£ = 
AiEXA' - A\ 2 A 2 + e 4p = 0, 



(35) 
(36) 



where a prime denotes differentiation with respect to p. 

If we introduce a new function D(p), where A(p) = din D/ dp, and choose A = 
— iE, Eq. ( |5ED simplifies to 

D" + -r^e 4p D = 0. (37) 



AE 2 



The general solution to this equation is given by 

D = Z 



1 2(a+0+) 
AE 



(3? 
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where Zq is an arbitrary, linear combination of ordinary, zero-order Bessel functions. 
Finally, Eq. fl35[) may now be solved by direct integration after substitution of Eq. 
■ We conclude, therefore, that one wavefunction satisfying Eq. QUID is 



^ (e 2 74£)l 1 exp 



(39) 



To summarize, we have considered the quantization of the homogeneous Bianchi 
class A cosmologies in the limit of high anisotropy. A number of exact solutions to 
the quantum constraints were found. We considered first the approach whereby the 
momentum constraints are reduced at the classical level prior to quantization. Wave- 
functions that represent tunnelling from a classically forbidden state were presented. 
These solutions may be expressed as a linear superposition of purely Euclidean wave- 
functions. We then considered a quantum reduction of the momentum constraints for 
the vacuum Bianchi IX cosmology. The modified Wheeler-DeWitt equation may be 
solved if appropriate conditions are satisfied. These solutions should prove useful for 
investigating some of the fundamental questions that arise in quantum cosmology. 

We would like to thank A. Higuchi for helpful communications on how the mo- 
mentum constraints can be satisfied at the quantum level. This work was supported 
by the Particle Physics and Astronomy Research Council (PPARC), UK. 
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